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Abstract 

The search for extraterrestrial intelligence (SETI) has been heavily influenced by 
solutions to the Drake Equation, which returns an integer value for the number of 
communicating civilisations resident in the Milky Way, and by the Fermi Paradox, 
glibly stated as: "If they are there, where are they?". Both rely on using average 
values of key parameters, such as the mean signal lifetime of a communicating 
civilisation. A more accurate answer must take into account the distribution of 
stellar, planetary and biological attributes in the galaxy, as well as the stochastic 
nature of evolution itself. This paper outlines a method of Monte Carlo realisation 
which does this, and hence allows an estimation of the distribution of key param- 
eters in SETI, as well as allowing a quantification of their errors (and the level of 
ignorance therein). Furthermore, it provides a means for competing theories of life 
and intelligence to be compared quantitatively. 

Keywords: Numerical, Monte Carlo, extraterrestrial intelligence, SETI, Drake 
Equation, Fermi Paradox 
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1 Introduction 

The science of SETI has always suffered from a lack of quantitative substance (purely 
resulting from its reliance on one-sample statistics) relative to its sister astronomical 
sciences. In 1961, Frank Drake took the first steps to quantifying the field by devel- 
oping the now-famous Drake Equation, a simple algebraic expression which provides 
an estimate for the number of communicating civilisations in the Milky Way. Unfor- 
tunately, its simplistic nature leaves it open to frequent re-expression, hence there are 
in fact many variants of the equation, and no clear canonical form. For the purpose of 
this paper, the following form will be used (Walters, Hoover & Kotra (1980) [27 1): 

N = R*f g f P n e fifif c L (1) 
With the symbols having the following meanings: 

N = The number of Galactic civilisations who can communicate with Earth 

i?* = The mean star formation rate of the Milky Way 

f g = The fraction of stars that could support habitable planets 

f p = The fraction of stars that host planetary systems 

n e = The number of planets in each system which are potentially habitable 

fl = The fraction of habitable planets where life originates and becomes complex 

fi = The fraction of life-bearing planets which bear intelligence 

f c = The fraction of intelligence bearing planets where technology can develop 

L = The mean lifetime of a technological civilisation within the detection window 

The equation itself does suffer from some key weaknesses: it relies strongly on mean 
estimations of variables such as the star formation rate; it is unable to incorporate the 
effects of the physico-chemical history of the galaxy, or the time-dependence of its 
terms. Indeed, it is criticised for its polarising effect on "contact optimists" and "con- 
tact pessimists", who ascribe very different values to the parameters, and return values 
of TV between icr 5 and 10 6 (!). 

A decade before, Enrico Fermi attempted to analyse the problem from a different angle, 
used order of magnitude estimates for the timescales required for an Earthlike civili- 
sation to arise and colonise the galaxy to arrive at the conclusion that the Milky Way 
should be teeming with intelligence, and that they should be seen all over the sky. This 
lead him to pose the Fermi Paradox, by asking, "Where are they?". The power of this 
question, along with the enormous chain of events required for intelligent observers to 
exist on Earth to pose it, has lead many to the conclusion that the conditions for life to 
flourish are rare, possibly even unique to Earth (Ward and Brownlee (2000) 11281 . The 
inference by Lineweaver (2001 ) [ 1 5 1 that the median age of terrestrial planets in the 
Milky Way is 1.8 ± 0.9 Gyr older than Earth would suggest that a significant number 
of Earthlike civilisations have had enough time to evolve, and hence be detectable: the 
absence of such detection lends weight to the so-called "rare Earth" hypothesis. How- 
ever, there have been many posited solutions to the Fermi Paradox that allow ETI to be 
prevalent, such as: 
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• They are already here, in hiding 

• Contact with Earth is forbidden for ethical reasons 

• They were here, but they are now extinct 

• They will be here, if Mankind can survive long enough 

Some of these answers are inherently sociological, and are difficult to model. Others 
are dependent on the evolution of the galaxy and its stars, and are much more straight- 
forward to verify. As a whole, astrobiologists are at a tremendous advantage in com- 
parison with Drake and Fermi: the development of astronomy over the last fifty years 
- in particular the discovery of the first extra solar planet (Mayor et al 1995 ifTTl ) and 
some hundreds thereafter, as well as the concepts of habitable zones, both stellar (Hart 
1979 irPD ) and galactic (Lineweaver et al 2004 lfl6l ) - have allowed a more in-depth 
analysis of the problem. However, the key issue still affecting SETI (and astrobiology 
as a whole) is that there is no consensus as to how to assign values to the key biological 
parameters involved in the Drake Equation and Fermi Paradox. Furthermore, there are 
no means of assigning confidence limits or errors to these parameters, and therefore no 
way of comparing hypotheses for life (e.g. panspermia - review, see Dose (1986)[9 | - 
or the "rare Earth" hypothesis (Ward and Brownlee (2000)[28|). This paper outlines 
a means for applying Monte Carlo Realisation techniques to investigate the parameter 
space of intelligent civilisations more rigorously, and to help assign errors to the result- 
ing distributions of life and intelligence. 

The paper is organised as follows: in section|2]the techniques are described; in section 
[3] the input data is discussed; in section H] the results from several tests are shown, and 
in section [5]the method is reviewed. 

2 Method 

The overall procedure can be neatly summarised as: 

1. Generate a galaxy of TV* stars, with parameters that share the same distribution 
as observations 

2. Generate planetary systems for these stars 

3. Assign life to some of the planets depending on their parameters (e.g. distance 
from the habitable zone) 

4. For each life-bearing planet, follow life's evolution into intelligence using stochas- 
tic equations 

This will produce one Monte Carlo Realisation (MCR) of the Milky Way in its en- 
tirety. The concept of using MCR techniques in astrobiology is itself not new: recent 
work by Vukotic & Cirkovic (2007,2008)1 24, 25 1 uses similar procedures to investigate 
timescale forcing by global regulation mechanisms such as those suggested by Annis 
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Figure 1: The stellar IMF used in this work (Scalo & Miller (1979) 122]). This is an 
example produced as part of a single MCR. 



(1999)[ D- In order to provide error estimates, this procedure must be repeated many 
times, so as to produce many MCRs, and to hence produce results with a well-defined 
sample mean and sample standard deviation. The procedure relies on generating pa- 
rameters in three categories: stellar, planetary and biological. 

2.1 Stellar Properties 

The study of stars in the Milky Way has been extensive, and their properties are well- 
constrained. Assuming the stars concerned are all main sequence objects allows much 
of their characteristics to be determined by their mass. Stellar masses are randomly 
sampled to reproduce the observed initial mass function (IMF) of the Milky Way, which 
is taken from Scalo & Miller (1979)022 (see FigureQ}- 

Stellar radii can then be calculated using (Prialnik 2000 |[T9l ): 



Rq 



Mr: 



(2) 



where n = 4 if the primary fusion mechanism is the p-p chain (M* < 1.1 Mq), and 
n = 16 if the primary fusion mechanism is the CNO cycle (M* > 1. 1M©). (Please 
note that in this paper, the subscript denotes the Sun, e.g. M Q indicates the value of 
one solar mass). The luminosity is calculated using a simple mass-luminosity relation: 



The main sequence lifetime therefore is: 



AT, 

Mr, 



(3) 
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Figure 2: The star formation history used in this work (Twarog (1980) H23V 
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The stars' effective temperature can be calculated, assuming a blackbody: 



TL = 



"I 1/4 



AirRlasB 



(4) 



(5) 



The star's age is sampled to reproduce the star formation history of the Milky Way 
(Twarog 1980 (23]), see Figureg] 

The metallicity of the star is dependent on the metallicity gradient in the galaxy, and 
hence its galactocentric radius, r ga i . This is sampled so that the surface mass density of 
the galaxy is equivalent to that of the Milky Way (assuming a simple two-dimensional 
disc structure): 



S(r ga/ ) = S e 



—r ga t/r h 



(6) 



where is the galactic scale length (taken to be 3.5 kpc). Therefore, given its galac- 
tocentric radius, the metallicity of the star, in terms of [^f\ , is calculated using the 
simple parametrisation of the abundance gradient: 



rad log 



r gai 

"gal,® 



(7) 



Where z gra d is taken to be 0.07 (Hou, Prantzos and Boissier (2000) [ 12|). The star is 
then placed at its galactocentric radius into one of four logarithmic spiral arms (mim- 
icking the four main logarithmic spiral arms of the Milky Way, see Figure|3|. 
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Figure 3: An example star map, taken from one MCR. Brighter areas indicate greater 
stellar mass 



Finally, stars are assigned planets based on their metallicity (sampling the current 
metallicity relations for extrasolar planets, see Figure [6}. Some stars will have no 
planets, some will have only one, and some may have multiple systems (based on the 
multiplicity of current extrasolar planetary systems 

2.2 Planetary Properties 

Since the discovery of 51 Peg B (Mayor et al 1995 [ 17]), the data garnered on exoplan- 
ets has grown at an increasing rate. At the time of writing, there are over 300 known ex- 
oplanets, discovered with a variety of observational techniques (Radial Velocity, Tran- 
sits, Microlensing, etc). This data provides distributions of planetary parameters that 
can be sampled from: the planetary mass function (PMF), the distribution of planetary 
orbital radii, and the host star metallicity distribution. 

Therefore, as with the stellar parameters, a population of planets can be created around 
the parent stars, with statistical properties matching what can be observed. However, 
this statistical data is still subject to strong observational bias, and the catalogues are 
still strongly incomplete. There is insufficient data to reproduce a distribution of terres- 
trial planets: therefore it is assumed that life evolves around the satellites of the planets 
simulated here. In essence, this constitutes a lower limit on the number of inhabited 
planets: the work of Ida and Lin (2004)[13| shows that, as a function of metallicity, 
habitable terrestrial planets are comparable in frequency (or higher) than currently de- 

1 Although not considered here, it should be noted that too high a metallicity may be as negative a factor 
as low metallicity, as high metallicity systems will tend to produce ocean worlds (Leger et al 2004 1141 1 
which may prove hostile to the formation of intelligent life. 
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Planet Mass (M'"') 

Figure 4: The planetary mass function, constructed from the Exoplanet Encyclopedia 
data (http://exoplanet.eu) 




Figure 5: The distribution of planetary orbital radii, constructed from the Exoplanet 
Encyclopedia data (http://exoplanet.eu) 
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Figure 6: The distribution of host star metallicity, constructed from the Exoplanet En- 



tectable giant planets. This data is hence still useful for illustrating the efficacy of the 
Monte Carlo method (at least, until observations of terrestrial exoplanets become sta- 
tistically viable). All this should be borne in mind when the results of this work are 
considered. 

Two further parameters to be sampled, for which observational data is difficult to ob- 
tain, are Puf e , the probability that life can originate on a given planet, and N resets , 
the number of biologically damaging events that a planet will experience during the 
evolution of its indigenous life (such as local supernovae, gamma ray bursts (Annis 
1999), cometary impacts, etc). At this point, the model goes beyond relatively well- 
constrained parameters, and becomes hypothesis. N resets is defined as a function of 
galactocentric radius - more resetting events (most notably SNe) occur as distance to 
the Galactic Centre decreases: 



Where /i re sets,o is set to 5, reflecting the so-called "Big Five" mass extinction events 
in the Earth's fossil record (Raup and Sepkoski (1982)). It should be noted that this 
is a simple parameterisation: the real distribution of resetting events is more complex, 
in particular for short GRBs and Type I SNe, which have less well-defined spatial 
relationships. This simple expression is provided to prevent a blurring of the issue: 
future work will include a more rigorous expression. All planetary data was collected 
from the Exoplanets Encyclopaedia website (http://exoplanet.eu): in total, 242 planets 
were used to construct the data, as these had entries in all required data fields. 



cyclopedia data (http://exoplanet.eu) 




(8) 
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2.3 Biological Parameters 
2.3.1 Life Parameters 

The model now enters the realm of essentially pure conjecture: all the available data 
for these parameters is derived from observations of a single biosphere, and hence there 
is little that can be done to constrain these parameters (at least without making wide- 
ranging assumptions about the mechanisms of life as a whole). The model implicitly 
assumes the "hard step" scenario of evolution (e.g. Carter 2007 [3|), i.e. life must 
achieve essential evolutionary goals in order to become intelligent organisms with the 
ability to construct sufficiently complex technological artifacts. The key biological 
parameters are listed as follows: 

• N stages'- The number of stages life must evolve through to become intelligent 

• Tj : The time required for each of these stages to be reached 

• Ti n t ■ The time required for an intelligent technological civilisation to form from 
life's creation 

• PannihUate '■ The probability that a reset event will cause complete annihilation 

The intelligence timescale Ti nt is calculated by the following stochastic process: if life 
does evolve on a planet, N sta ges is randomly sampled. The resetting events N reset8 
are placed uniformly throughout each of the stages. If N resets > N stag es, then any 
given stage may suffer several reset events. Then the following procedure occurs for 
each stage i; 

1 . Tj is sampled 

2. If a reset occurs, test if that reset results in annihilation: if annihilation occurs, 
life is exterminated, and the process ends; otherwise, i decreases by 1. 

3. Ti n t is increased by some fraction of r, (or simply by Tj if no resets occur) 

4. Increase i by 1, and return to 1 

This procedure continues until either a) life has reached the final stage (an intelligent 
technologically capable civilisation has evolved), or b) the intelligence timescale be- 
comes greater than the main sequence lifetime of the parent star: Ti nt > tms- Once a 
civilisation has formed, it is assumed that detectable signals or signal leakage begins to 
be emitted. This emission will continue until either (a) the civilisation destroys itself 
or (b) the parent star moves off the main sequence (see next section). This in itself is 
a conservative estimate of the length of time a civilisation may be detected, as civilisa- 
tions may in fact migrate between the stars, or use stellar engineering to prolong their 
parent's stars life. 

The P an nihilate parameter is defined as: 



P .... . = 1 _ p r ga i-r ga i tQ 
1 annihilate — ± c 



(9) 
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Along with the metallicity gradient and the N resets parameter specified in the pre- 
vious sections, this defines a Galactic Habitable Zone which mimics that defined by 
Lineweaver et al (2004)|[T5|: an annular region of the galaxy where the metallicity 
is sufficiently high for planets (and life) to form, and the number of biologically de- 
structive events is sufficiently small (and their destructive capability sufficiently low) 
to allow life to evolve at all. 



2.3.2 Civilisation Parameters 

Once a technologically capable civilisation has formed, it must move through a "fledgling 
phase": it is susceptible to some catastrophic event caused partially or fully by its 
own actions (e.g. war, plague, catastrophic climate change, botched macro-engineering 
projects). This is described by the parameter Pdestroy- the probability that a fledgling 
civilisation will destroy itself. If a civilisation can survive this phase, it becomes suf- 
ficiently advanced to prevent such self-destruction events, and becomes stable, on a 
timescale T a d v - If a civilisation is destroyed, then it will survive some fraction of r a dv 
before destruction. 

What advanced civilisations can then do is at the behest of the user: civilisations may 
colonise all planets within their solar system, resulting in signals appearing on all plan- 
ets in that system. Probes may be sent into the galaxy at large, which could define an ex- 
plorable volume of the galaxy for a given advanced civilisation. Civilisations may even 
attempt to generate new biospheres on neighbouring planets - the "directed pansper- 
mia" model (Crick and Orgel 1972 [7 1). The results in this paper assume colonisation 
of all uninhabited planets in the system will occur once a civilisation becomes ad- 
vanced (colonisation being an umbrella term for placing both manned and unmanned 
installations on the surface or in orbit of said planets 

The key civilisation parameters are: 

• T a d v - The timescale for a civilisation to move from "fledgling" to "advanced". 

• Pdestroy- The probability that a fledgling civilisation will destroy itself. 

• L S ig n ai : The lifetime of any signal or leakage from a civilisation 
The signal lifetime of a self-destroying civilisation is: 

L signal — -ET~ a dv (10) 

Where a; is a uniformly sampled number between and 1 . If the civilisation becomes 
advanced, this becomes 

2 It is reasonable to assume that a pre-existing biosphere, perhaps with intelligent life would be a more 
suitable candidate for colonisation rather than a lifeless planet: however, to prevent statistical confusion (and 
to avoid dealing with the possibilities of interplanetary conflict between intelligent species) colonisation 
occurs only on uninhabited planets 
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L signal — T~MS ~ T int (11) 

(i.e. civilisations exist until their parent star leaves the Main Sequence). For a planet 
colonised by an advanced civilisation, this is 

L signal = T MS ~ T int ~ T adv (12) 

At the end of any MCR run, each planet will have been assigned a habitation index 
based on its biological history. 



^-inhabit 



— 1 Biosphere which has been annihilated 

Planet is lifeless 

1 Planet has primitive life 

2 Planet has intelligent life (13) 

3 Planet had intelligent life, but it destroyed itself 

4 Planet has an advanced civilisation 

5 Planet has been colonised by an advanced civilisation 



If the planet has habitation index 1 or 2, the biological process has been ended by the 
destruction of the parent star. Planets with an index of -1 or 1 may contain biomarkers 
in their atmosphere (e.g. ozone or water spectral features) which could be detected. 
Planets with an index of 2 or higher will emit signals or signal leakage. Planets with 
an index of 4 or 5 may display evidence of a postbiological civilisation or of large 
scale "macro-engineering" projects , e.g. Dyson spheres (Dyson 1960[10|). Signals 
from these systems may even be consistent with those expected from Kardashev Type 
II civilisations (those which can harness all the energy of their parent star), and hence 
could produce characteristic stellar spectral signatures that Earth astronomers could 
detect. 



3 Inputs 

When comparing hypotheses, it is important to keep the number of free parameters to 
a minimum. To achieve this, a biological version of the Copernican Principle is in- 
voked: the Terran biosphere is not special or unique. Therefore, life on other planets 
will share similar values of characteristic parameters. This implies that many parame- 
ters can be held constant throughout all hypotheses: in particular, those that deal with 
the approach to intelligence. These are sampled from Gaussian distribution functions 
(assuming that each variable is the result of many subfactors, and implicitly applying 
the Central Limit Theorem). The mean value for N sta ges is taken to be 6, reflecting the 
major stages life went through on Earth - biogenesis, the advent of bacteria, the advent 
of eukaryotes, combigenesis, the advent of metazoans, and the birth of technological 
civilisation (Carter 2007 O). 100 MCRs were produced for each hypothesis: the pa- 
rameters used for all testing can be found in Table [J The three hypotheses tested are 
described below. 
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3.1 The Panspermia Hypothesis 

This well-documented theory suggests life may spread from one originating planet to 
many others, causing life to form concurrently in multiple systems. In this model, if 
life forms on any planet in a star system, then other planets may be seeded, according 
to the following prescription: 



Life can move between planets if the distance between the destination planet and the 
original inhabited planet (Ai?) is sufficiently small, and the stellar luminosity is suffi- 
ciently low that ionising radiation will not destroy biological organisms in transit. It is 
assumed that life emerges on the original planet while a significant amount of planetary 
bombardment is in progress, and the seeding of other planets (from rocky fragments ex- 
pelled from the inhabited planet's surface by asteroid impacts) is therefore effectively 
instantaneous. Pdestroy is specified a priori to be 0.5, reflecting current ignorance 
about the development of extraterrestrial civilisation. 

This hypothesis restricts panspermia to planets in the same star system only: unfortu- 
nately, the model is currently incapable of modelling the interesting possibility of in- 
terstellar panspermia (Napier 2004 [18], Wallis and Wickramasinghe 2004 [26|). This 
in itself is strong motivation to improve the model further, and is left as an avenue for 
future work. 



Due to the current dearth of data on Earth-like exoplanets, the famous "Rare Earth 
Hypothesis" (Ward and Brownlee (2000)[28|) cannot be tested fully. A different hy- 
pothesis is possible (assuming Jovian planets in the habitable zone may have stable 
Earth-like moons on which life could arise) with this code. The hypothesis is as fol- 
lows: once life arises, it is assumed to be tenacious, and can easily evolve into intel- 
ligence: however, the initial appearance of life is itself hard: the planetary niche must 
satisfy certain criteria before life can evolve. The criteria are: 

1 . The planet is in the stellar habitable zone 

2. The parent star mass is less than 2M Q 

3. The system has at least one extra planet (protection from bombardment) 

4. The star's metallicity must be at least solar or higher (Z* > Zq) 

If all criteria are met, then life can evolve. Pdestroy is set to 0.5 for the same reason 
as previously. It should be noted again that this is significantly different from the Rare 
Earth hypothesis, and should not be considered as a comparison. A more accurate 
comparison would have the last few stages of life (the development of metazoans and 




if planet is in stellar habitable zone 

if any planet in the same solar system is inhabited 



(14) 



3.2 The Rare Life Hypothesis 



A Testbed for Hypotheses ofETI 



13 



intelligent life) be rare: this would be better suited for future work after the model has 
been refined. 

3.3 The Tortoise and Hare Hypothesis 

As a final example, this hypothesis demonstrates the ability of the code to model more 
sociological hypotheses. Life evolves easily on many planets, but the evolution towards 
intelligence and advanced civilisation is more complex. Essentially, civilisations which 
arise too quickly are more susceptible to self-destruction, whereas those which take 
longer to emerge are more likely to survive the fledgling phase. This is parametrised 
as: 

Pdestroy — (15) 
Tint 

Where to = N s tagesT m i n is a normalisation constant which describes the minimum 
evolution timescale (if N resets = 0, and each stage takes the minimum amount of time 
Tmin) (Note:If this calculation is applied to Mankind, our own destruction probability 
is around 0.8 (!)). If a planet is in the habitable zone, then Puf e = 1, in the same 
vein as the other hypotheses. It should be noted that this is an oversimplification: the 
potential existence of liquid water on Enceladus and Europa (both outside the solar 
habitable zone) suggests that Pu f e should be a distribution in orbital radius with tails 
that extend outside the zone: however, for the purposes of this paper, a simple step 
function will be sufficient. 



Table 1 : Parameters used for all hypotheses 



Hypothesis 


Parameter 


Mean 


Standard Deviation 


All 


A stages 


6 


1 


All 


n 


0.8 Gyr 


0.25 Gyr 


All 


'adv 


2.5 xlO" 4 Gyr 


1.0 xl0" 4 Gyr 


Hypothesis 


Parameter 


Value 




All 


Nmcr 


100 




Panspermia 


Plife 


as equation d 1 4b 




Panspermia 


* destroy 


0.5 




Rare Life 


Plife 


1 (if all criteria met) 




Rare Life 


^destroy 


0.5 




Tortoise & Hare 


Plife 


1 (if in stellar habitable zone) 




Tortoise & Hare 


Pdestroy 


as equation (TT5b 
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4 Results 

4.1 The Panspermia Hypothesis 

The single number statistics for this hypothesis can be seen in Table |2] As expected, 
around half of all emerging civilisations destroy themselves (as Pdestroy = 0.5). Around 
0. 1 % of all planets are inhabited by either primitive or intelligent life during the Galaxy's 
lifetime, and roughly 10% of all inhabited planets host intelligent lifeforms. Compared 
to the Tortoise and Hare hypothesis (see below), the inhabited fraction has not greatly 
increased with possibility of seeding amongst neighbours: this is mostly due to the low 
multiplicity of current exoplanets (around 11%). 



Table 2: Statistics for the Panspermia Hypothesis 

Variable Mean Standard Deviation 

anets,total 

4.7295 x 10 s 401530 

NrnhaUted 690983.63 53 

Nfiedgiing 75923.04 8 

Ndestroyed 37958.07 11 

Nadvanced 37964.97 20 



Figure |7] shows the Galactic Habitable Zone for this hypothesis: at low galactocentric 
radii, the annihilation rate increases to the point where civilisations cannot develop (al- 
though life itself can emerge there, it is quickly eradicated by destructive astrophysical 
events). The right panel of Figure|7]shows the breakdown of inhabited planets by type 
over the life of the Galaxy. The large number of planets with index -1 (annihilated bio- 
spheres) indicates the solar systems with low galactocentric radii which are metal rich, 
and form planets (and hence life) more easily, but suffer rapid biological catastrophe. 
The colonisation fraction (planets with habitation index 5) is quite low: this is most 
likely due to panspermia seeding reducing the number of neighbouring empty planets 
available for colonisation. 

The distribution of stellar properties for civilised planets (Figure|8]l indicates a predilec- 
tion for low mass stars. This comes from the selection bias of the so-called "Hot 
Jupiters": these objects are of such low orbital radii that if they are to exist in the hab- 
itable zone, it must also be of low radius, and hence belong to a low-mass star. The 
metallicity reflects the Galactic Habitable Zone: civilisations have a minimum galac- 
tocentric radius, and hence a maximum metallicity. Considering the signal lifetime 
distribution, it is heavily influenced by the main sequence lifetime of the parent stars, 
which is a function of mass. Studying the mass relation, it is clear that the two are cor- 
related. If the signals are plotted as a function of time (Figure|9] right panel), it can be 
seen that the signal history undergoes a transition from low to high N around t ~ tn, 
where tn is the Hubble time (i.e., the age of the universe). This is symptomatic of the 
Copernican Principle invoked to constrain the majority of the biological parameters: 
on average, any given civilisation will have a intelligence timescale similar to ours. 
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r .-' ( k P c ) Habitation Index 

Figure 7: The distribution of galactocentric radius (left) and habitation index (right) 
under the Panspermia Hypothesis. 




M. (M„) Z. [Fe/H] 

Figure 8: The distribution of host star mass (left) and metallicity (right) under the 
Panspermia Hypothesis. 




Figure 9: The distribution of signal lifetimes (left) and the signal history of the Galaxy 
( right) under the Panspermia Hypothesis. 
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Habitation nde 



Figure 10: The distribution of galactocentric radius (left) and habitation index (right) 
under the Rare Life Hypothesis. 



4.2 The Rare Life Hypothesis 

The stringent nature of this hypothesis can be seen in its decreased population of life- 
bearing planets in comparison to the other hypotheses (with an inhabited fraction of 
only 0.01%). However, despite this stringency, there is a significant number of civilisa- 
tions forming over the Galaxy's lifetime: around 1 % of all life-bearing planets produce 
intelligent species. Again, as expected around half of all civilisations destroy them- 
selves. 



Table 3: Statistics for the Rare Life Hypothesis 

Variable Mean Standard Deviation 

N planets ,total 4.770 x 10* 381 

habited 80090 22 

Nfledgling 728.6 1 

Ndestroyed 367.3 

N advanced 361.2 2 



The distribution of galactocentric radii (Figure[T0l left panel) is quite sharply peaked: 
this is because life requires Z* > Zq (which is clearly seen in FigureQl] right panel) 
under this hypothesis, and fixes a maximum r ga i , beyond which life cannot begin. Con- 
sidering the habitation index (FigureQl)] right panel), it can be seen that there is a great 
deal of colonisation occurring. This is due to the criterion that life can only develop in 
multiple planetary systems, and hence any intelligence which becomes advanced must 
take part in colonisation. The distribution of star mass, although greatly reduced, is 
qualitatively the same as for the Panspermia Hypothesis, and hence the signal lifetime 
has a similar distribution also. 
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Figure 12: The distribution of signal lifetimes (left) and the signal history of the Galaxy 
( right) under the Rare Life Hypothesis. 
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Figure 13: The distribution of galactocentric radius (left) and habitation index (right) 
under the Tortoise and Hare Hypothesis. 

4.3 The Tortoise and Hare Hypothesis 

The biological data for this hypothesis is not markedly different from the Panspermia 
Hypothesis, again with an inhabited fraction of around 0.1%. However, the true differ- 
ence is in the civilisation data: instead of there being an approximate 1 : 1 ratio between 
self-destruction and advancement, it can be seen that this hypothesis weakly favours 
self-destruction for a given fledgling civilisation, with a ratio of approximately 1.38:1. 
Despite this, around 4% of inhabited planets develop intelligent life that survives the 
fledgling phase. 



Table 4: Statistics for the Tortoise and Hare Hypothesis 



Variable 


Mean 


Standard Deviation 


N planets , total 


4.770 x 10 s 


1593 


N inhabited 


684399.26 


2 


N fledgling 


75200.3 


20 


N destroyed 


43626.82 


1 


^advanced 


31573.52 


20 



The Galactic Habitable Zone is similar to that for the Panspermia Hypothesis: the weak 
favouring of self-destruction can be seen in the right panel of Figure[l3] Comparing the 
data between this hypothesis and the Panspermia Hypothesis, both the qualitative and 
the quantitative behaviour is similar. This reflects the similarity of the two hypotheses 
in terms of biological parameters, and that the key difference is in the sociological 
parameter P d estroy 
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Figure 14: 77ze distribution of host star mass (left) and metallicity (right) under the 
Tortoise and Hare Hypothesis. 




Figure 15: The distribution of signal lifetimes (left) and the signal history of the Galaxy 
( right) under the Tortoise and Hare Hypothesis. 
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5 Conclusions 

The results for all three hypotheses show clear trends throughout. In all cases, the 
inhabited planets orbit low mass stars (a symptom of the Hot Jupiter bias present in 
current data); the signal lifetimes are correspondingly dependent on these masses (as 
expected, given its functional form). The signal history of all three hypotheses shows 
a period of transition between low signal number and high signal number at around 
t ~ tu, again symptomatic of the Biological Copernican Principle used to constrain 

7"is ^stages-) etc. 

This paper has outlined a means by which key SETI variables can be estimated, taking 
into account the diverse planetary niches that are known to exist in the Milky Way and 
the stochastic evolutionary nature of life, as well as providing estimates of errors on 
these variables. However, two notes of caution must be offered: 

1. The reader may be suspicious of the high precision of the statistics quoted: it is 
worth noting that the standard deviations of these results are indeed low, and the 
data is precise: but, its accuracy is not as certain. The output data will only be as 
useful as the input data will allow (the perennial "garbage in, garbage out" prob- 
lem). Current data on exoplanets, while improving daily, is still insufficient to 
explore the parameter space in mass and orbital radii, and as such all results here 
are very much incomplete. Conversely, as observations improve and catalogues 
attain higher completeness, the efficacy of the Monte Carlo Realisation method 
improves also. Future studies will also consider planetary parameters which are 
sampled as to match current planet formation theory, rather than current obser- 
vations. 

2. The method currently does not produce a realistic age metallicity relation (AMR). 
Age, metallicity and galactocentric radius are intrinsically linked in this setup: 
to obtain realistic data for all three self-consistently requires an improved three- 
dimensional galaxy model which takes into account its various components (the 
bulge, the bar, etc), as well as the time evolution of the Galaxy. Future work will 
attempt to incorporate a more holistic model which allows all three parameters to 
be sampled correctly. In particular, future efforts will be able to take advantage 
of better numerical models for the star formation history (e.g. Rocha-Pinto et al 
2000 [21 1), and the spatial distribution of stars (e.g. Dehnen and Binney 1998 

El). 

Although this paper applies the "hard step scenario" to the biological processes mod- 
elled, the method itself is flexible enough to allow other means of evolving life and 
intelligence. The minutiae of how exactly the biological parameters are calculated 
do not affect the overall concept: this work has shown that it is possible to simulate 
a realistic backdrop (in terms of stars and planets) for the evolution of ETI, whether 
it is modelled by the "hard step" scenario or by some other stochastic method. In- 
corporating new empirical input from the next generation of terrestrial planet finders, 
e.g. Kepler (Borucki et al 2008 [2]), as well as other astrobiological research into the 
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model, alongside new theoretical input by adding more realistic physics and biology 
will strengthen the efficacy of this Monte Carlo technique, providing a new avenue of 
SETI research, and a means to bring many disparate areas of astronomical research 
together. 
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